Available  online  at  www.sciencedirect.com 


ELSEVIER 


*V  ScienceDirect 


Journal  of  Power  Sources  162  (2006)  388-399 


JOURNAL  Of 


www.elsevier.com /locate /jpowsour 


Model-based  condition  monitoring  of  PEM  fuel  cell  using 

Hotelling  T 2  control  limit 

X.  Xuea,  J.  Tang  a’ * ,  N.  Sammesa,  Y.  Dingb 

a  Department of  Mechanical  Engineering,  University  of  Connecticut,  191  Auditorium  Road,  Unit  3139,  Storrs,  CT  06269-3139,  United  States 
b  Department  of  Industrial  Engineering,  Texas  A&M  University,  College  Station,  IX  77843-3131,  United  States 

Received  30  May  2006;  received  in  revised  form  30  June  2006;  accepted  3  July  2006 

Available  online  23  August  2006 


Abstract 

Although  a  variety  of  design  and  control  strategies  have  been  proposed  to  improve  the  performance  of  polymer  electrolyte  membrane  (PEM)  fuel 
cell  systems,  temporary  faults  in  such  systems  still  might  occur  during  operations  due  to  the  complexity  of  the  physical  process  and  the  functional 
limitations  of  some  components.  The  development  of  an  effective  condition  monitoring  system  that  can  detect  these  faults  in  a  timely  manner  is 
complicated  by  the  operating  condition  variation,  the  significant  variability/uncertainty  of  the  fuel  cell  system,  and  the  measurement  noise.  In  this 
research,  we  propose  a  model-based  condition  monitoring  scheme  that  employs  the  Hotelling  T 2  statistical  analysis  for  fault  detection  of  PEM  fuel 
cells.  Under  a  given  operating  condition,  the  instantaneous  load  current,  the  temperature  and  fuel/gas  source  pressures  of  the  fuel  cell  are  measured. 
These  measurements  are  then  fed  into  a  lumped  parameter  dynamic  fuel  cell  model  for  the  establishment  of  the  baseline  under  the  same  operating 
condition  for  comparison.  The  fuel  cell  operation  is  simulated  under  statistical  sampling  of  parametric  uncertainties  with  specified  statistics  (mean 
and  variance)  that  account  for  the  system  variability/uncertainty  and  measurement  noise.  This  yields  a  group  of  output  voltages  (under  the  same 
operating  condition  but  with  uncertainties)  as  the  baseline.  Fault  detection  is  facilitated  by  comparing  the  real-time  measurement  of  the  fuel  cell 
output  voltage  with  the  baseline  voltages  by  employing  the  Hotelling  T 2  statistical  analysis.  The  baseline  voltages  are  used  to  evaluate  the  output 
T2  statistics  under  normal  operating  condition.  Then,  with  a  given  confidence  level  the  upper  control  limit  can  be  specified.  Fault  condition  will 
be  declared  if  the  T2  statistics  of  real-time  voltage  measurement  exceeds  the  upper  control  limit.  This  model-based  robust  condition  monitoring 
scheme  can  deal  with  the  operating  condition  variation,  various  uncertainties  in  a  fuel  cell  system,  and  measurement  noise.  Our  analysis  indicates 
that  this  scheme  has  very  high  detection  sensitivity  and  can  detect  the  fault  conditions  at  the  early  stage. 

©  2006  Elsevier  B.Y.  All  rights  reserved. 
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1.  Introduction 

Polymer  electrolyte  membrane  (PEM)  fuel  cell  systems  oper¬ 
ate  at  relatively  low  temperatures  and  are  promising  candidates 
for  future  clean  energy  sources.  These  systems  are  intrinsically 
complex,  with  coupling  effects  of  fluid,  heat,  electrochemical 
reaction,  phase  change,  etc.  There  are  also  physical  property 
limitations  for  fuel  cell  components  such  as  membrane  and  elec¬ 
trodes.  Therefore,  fault  conditions  and  even  failure  may  occur 
under  practical  operating  conditions.  In  fact,  the  reliability  and 
durability  are  considered  the  main  challenges  to  the  fuel  cell 
research  community. 
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The  fault  conditions,  which  might  be  temporary  or  perma¬ 
nent,  usually  are  closely  related  to  a  critical  component  of  the 
fuel  cell,  the  membrane-electrode  assembly  (ME A).  Typically, 
three  fault  modes  may  be  involved  in  fuel  cell  failures,  i.e.,  the 
dehydration  and  drying  of  the  membrane,  the  fuel/gas  starvation 
of  electrochemical  reaction,  and  the  leak  of  the  membrane.  It  is 
well  known  that  the  electrolyte  membrane  needs  to  be  appro¬ 
priately  hydrated  in  order  to  efficiently  conduct  the  hydrogen 
protons  and  prevent  the  occurrence  of  localized  hot  spots.  The 
dehydration  and  drying  of  the  membrane  lead  to  the  increase  of 
the  internal  resistance  and  a  larger  output  voltage  loss,  which  in 
turn  raises  the  local  temperature  of  the  membrane.  If  the  local 
overheating  lasts  for  a  long  period  of  time,  hot  spot  may  occur 
and  eventually  permanently  damage  the  membrane.  Water  trans¬ 
port  across  the  membrane  and  water  content  therein  depend  on 
the  combinational  effects  of  the  electro-osmotic  drag  force  on 
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Nomenclature 

a 

water  activity 

Amem 

effective  area  of  membrane  (m2) 

Ci 

specific  heat  (J/(kg  K)-1) 

C 

water  concentration  (molm-3) 

Av 

water  diffusion  coefficient  (m2  s-1) 

F 

Faraday’s  constant  (A  s  mol)-1 

^latent 

latent  heat  of  water  phase  change  (J  kg-1) 

M 

effective  convection  heat  transfer  coefficient 

i 

current  density  (A  m-2) 

I 

load  current  (A) 

J w 

mole  flux  of  water  transport  in  the  membrane 
(mol(m2  s)-1) 

Ki 

flow  rate  coefficient  (kg(Pas)-1) 

m 

species  mass  in  the  channels  (kg) 

M 

species  mole  mass  (kg  mol-1) 

4Tm,dry 

mole  mass  of  dry  membrane  (kg  mol-1) 

nd 

electro-osmotic  drag  coefficient 

N 

number  of  cells  in  the  stack 

P 

species  partial  pressure  (Pa) 

fin 

membrane  thickness  (m) 

T 

temperature  (K) 

V 

channel  volume  (m3) 

W 

flow  rate  (kg  s-1) 

Greek  symbols 

p 

humidity  ratio 

0 

relative  humidity 

A 

water  content 

Pm,  dry 

dry  density  of  membrane  (kgm-3) 

r 

time  constant  (s) 

& 

empirical  coefficients  of  activation  overvoltage 

Superscripts  and  subscripts 

an 

anode 

atm 

ambient  air 

body 

cell  body 

ca 

cathode 

cw 

cooling  water 

hum 

humidification 

in 

inlet 

mem 

membrane 

02 

oxygen 

out 

outlet 

phase 

phase  change 

s 

source 

sat 

saturation 

V 

vapor 

w 

water 

the  water  molecules  and  water  back-diffusion  from  the  cathode 
side  to  the  anode  side.  It  is  believed  that  the  electro-osmosis  drag 
effect  is  larger  than  that  of  the  back-diffusion,  and  the  water  con¬ 
tent  of  the  membrane  is  mainly  dependent  on  that  of  the  anode 
side  [1-3].  If  the  fuel  is  not  adequately  humidified,  the  mem¬ 


brane  may  be  dehydrated  at  the  anode  side,  while  flooding  may 
occur  at  the  cathode  side  due  to  water  generated  by  the  electro¬ 
chemical  reaction.  The  fuel/gas  starvation  might  have  several 
possible  causes  [1-4].  The  first  stage  may  be  caused  by  the 
channel  flow  variation  within  the  stack  [4],  i.e.,  the  flow  resis¬ 
tance  directly  resulting  in  the  fuel/gas  starvation  in  the  channels, 
and  consequently  the  electrolyte  fuel/gas  starvation.  The  reasons 
include  the  liquid  water  droplets  forming  in  the  flow  channels, 
the  temperature  variation,  and  geometry  deviation.  The  second 
stage  fuel/gas  starvation  may  be  caused  by  the  electrode  pores 
blocked  by  the  liquid  water,  which  is  termed  “flooding”.  This 
phenomenon  generally  occurs  at  the  cathode  side.  The  fuel/gas 
starvation  may  interrupt  the  electrochemical  reaction  and  lead  to 
a  rapid  loss  in  the  output  voltage.  In  the  worst  situation,  it  may 
even  lead  to  the  decomposition  of  the  fuel  cell  component  and 
permanently  damage  the  cells  [3,4].  Another  fault  mode  is  the 
leak  of  the  membrane,  which  is  due  to  the  fracture  and/or  hole 
of  the  membrane  [4] .  While  the  holes  in  the  membrane  may  be 
caused  by  the  hot  spot,  the  fracture  of  the  membrane  is  directly 
resulted  from  the  mechanical  stress  concentration.  Under  the 
dynamic  operating  conditions,  the  pressure  difference  across  the 
ME  A  may  break  the  membrane  [4] . 

In  order  to  improve  the  reliability  and  the  overall  perfor¬ 
mance  of  PEM  fuel  cell  systems,  a  number  of  design  and  control 
strategies  have  been  proposed  and  some  examples  can  be  found 
in  Refs.  [5-8].  Currently,  however,  there  exists  very  limited 
research  on  the  monitoring  and  online  fault  detection  of  PEM 
fuel  cells.  Recently,  Hissel  et  al.  [9]  explored  a  fuzzy  diagnosis 
method  for  PEM  fuel  cells.  A  Sugeno-type  fuzzy  model  of  a 
PEM  fuel  cell  was  developed.  The  fuzzy  model  and  the  phys¬ 
ical  PEM  fuel  cell  process,  subjected  to  the  same  input,  were 
placed  in  parallel  to  each  other,  and  the  output  voltage  differ¬ 
ence  between  the  model  predication  and  the  measurement  was 
used  to  infer  the  fault  occurrence.  They  considered  the  cases 
of  accumulation  of  nitrogen  and  water  in  the  anode  compart¬ 
ment  and  the  drying  of  the  membrane.  The  main  difficulty  in 
the  online  monitoring  for  PEM  fuel  cells  stems  from  the  sys¬ 
tem  complexity.  While  over  the  years  PEM  fuel  cells  have  been 
modeled  at  various  levels  with  different  focuses,  a  global  model 
capable  of  characterizing  the  dynamic  and  transient  behavior  of 
the  fuel  cells  is  still  being  pursued.  Moreover,  the  modeling  of 
certain  components  is  based  on  phenomenological  characteriza¬ 
tions  and  the  fuel  cell  system  has  inherent  variability /uncertainty 
due  to  manufacturing  tolerance,  etc.  These  issues,  together  with 
the  online  measurement  noise,  make  it  very  difficult  to  develop 
a  deterministic  model-based  monitoring  system.  For  example, 
a  direct  level-difference  comparison  between  the  online  mea¬ 
surement  and  a  deterministic  baseline  predicted  by  the  dynamic 
model  of  a  healthy  fuel  cell  might  lead  to  significant  number  of 
false  alarms.  Meanwhile,  a  fuel  cell  system  normally  operates 
under  dynamic  and  varying  conditions,  e.g.,  constant  changes  in 
load  current,  temperature,  flow  rate,  etc.  Therefore,  one  would 
need  to  collect  a  huge  amount  of  baseline  data  under  all  these 
varying  conditions  if  a  non-model  based  scheme  is  used,  which 
would  be  extremely  costly  and  even  infeasible. 

Building  upon  the  previous  studies  and  understandings  men¬ 
tioned  above,  in  this  paper  we  explore  the  development  of  a 
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model-based  health  monitoring  strategy  for  PEM  fuel  cells  using 
the  statistical  analysis.  This  proposed  strategy  can  deal  with  the 
operating  condition  change  and  system  uncertainty /noise  simul¬ 
taneously.  Under  a  given  operating  condition,  the  instantaneous 
load  current,  the  temperature  and  fuel/gas  source  pressures  of 
the  PEM  fuel  cell  are  measured.  These  measurements  are  then 
fed  into  a  lumped  parameter  dynamic  fuel  cell  model  for  the 
establishment  of  the  baseline  under  the  same  operating  condi¬ 
tion  for  comparison.  In  order  to  account  for  the  system  vari¬ 
ability/uncertainty,  the  fuel  cell  operation  is  simulated  under 
statistical  sampling  of  parametric  uncertainties  with  specified 
statistics  (mean  and  variance).  This  yields  a  group  of  output  volt¬ 
ages  (under  the  same  operating  condition  but  with  uncertainties) 
as  the  baseline.  Fault  detection  is  facilitated  by  comparing  the 
real-time  measurement  of  the  fuel  cell  output  voltage  with  the 
baseline  voltages  from  the  statistical  analysis  standpoint.  This 
proposed  strategy  involves  two  aspects  of  advancement.  Build¬ 
ing  upon  the  recent  progress,  we  first  establish  a  system-level 
dynamic  model  of  PEM  fuel  cell  to  characterize  the  compli¬ 
cated  interactions  of  the  temperature,  gas  flow,  phase  change  in 
the  anode  and  cathode  channels,  and  membrane  humidification 
under  operating  conditions.  In  this  model,  we  explicitly  take 
into  account  the  phase  change  effect  so  that  the  water  content 
and  heat  inside  the  PEM  fuel  cell,  which  are  critical  to  the  nor¬ 
mal  operation  of  the  fuel  cell,  can  be  properly  evaluated.  At 
the  fault  detection  stage,  we  incorporate  the  Hotelling  T 2  statis¬ 
tics  method  [10]  into  the  condition  monitoring.  The  Hotelling 
T2  statistics  is  a  powerful  tool  of  multi- variate  statistical  anal¬ 
ysis,  and  has  been  widely  used  in  process  and  quality  controls 
[11,12].  Essentially,  the  T 2  value  of  the  response  data  is  used  to 
measure  the  overall  conformance  to  an  established  standard.  In 
the  proposed  condition  monitoring  scheme  for  PEM  fuel  cell,  the 
aforementioned  baseline  voltages  obtained  from  the  simulations 
under  the  same  input  are  used  to  evaluate  the  output  T2  statistics 
under  normal  operating  condition.  Then,  with  a  given  confidence 
level  the  upper  control  limit  can  be  specified.  Fault  condition  will 
be  declared  if  the  T2  statistics  of  real-time  voltage  measurement 
exceeds  the  upper  control  limit.  A  series  of  numerical  studies  are 
carried  out  to  demonstrate  the  system  performance.  Our  analysis 
shows  that  this  proposed  health  monitoring  system  is  capable  of 


detecting  fuel  cell  damage  under  complicated  environment  with 
system  variability/uncertainty  and  measurement  noise. 

2.  Health  monitoring  system  development 

2.7.  Condition  monitoring  structure 

The  structure  of  the  proposed  health  monitoring  system  is 
illustrated  in  Fig.  1 .  The  main  idea  is  to  use  the  statistical  analysis 
to  compare  the  online  measurement  of  the  fuel  cell  output  volt¬ 
age  to  the  baseline  predicted  by  the  healthy  reference  model.  The 
reference  model  is  placed  in  parallel  to  the  physical  fuel  cell  and 
is  subjected  to  the  same  inputs.  More  specifically,  the  physical 
fuel  cell  variables  such  as  fuel/gas  source  pressures,  environment 
temperature,  and  load  current  will  be  measured  using  various 
sensors.  All  these  measurement  data  will  be  used  as  inputs  to 
the  reference  model.  The  healthy  fuel  cell  output  voltage  predic¬ 
tion  under  the  specific  operating  condition  is  then  obtained  using 
the  fuel  cell  mathematical  model.  An  important  feature  of  the 
present  research  is  that  the  system  uncertainties/measurement 
noise  will  be  taken  into  account.  Generally,  the  system  uncer¬ 
tainties  stem  from  the  modeling  error  as  well  as  the  cell-to-cell 
difference  due  to  manufacturing  tolerance.  Mathematically,  the 
system  uncertainties  can  be  characterized  by  perturbations  (with 
certain  mean  and  variance)  to  physical  parameters  involved  in 
the  fuel  cell  model.  Therefore,  in  simulating  the  healthy  fuel 
cell  output  voltage,  we  will  run  the  reference  model  multi¬ 
ple  times  with  statistical  sampling  of  parametric  uncertainties. 
A  group  of  output  voltages  (with  uncertainties/noise)  can  be 
obtained  as  the  healthy  fuel  cell  baseline.  Fault  detection  is  then 
facilitated  by  comparing  the  real-time  measurement  of  the  fuel 
cell  output  voltage  with  the  baseline  voltages  by  employing  the 
Hotelling  T2  statistical  analysis.  The  Hotelling  analysis  includes 
two  stages,  the  baseline  construction  and  the  statistical  analysis 
evaluation  of  the  online  measurement  (monitored)  data.  At  the 
phase  one  stage,  the  voltage  outputs  of  the  healthy  fuel  cell  pre¬ 
dicted  by  multiple  simulations  with  system  uncertainties/noise 
will  be  analyzed  to  establish/confirm  the  control  limit  under  a 
given  confidence  level.  A  self-checking  procedure  is  employed 
to  either  adjust  the  control  limit  (and  the  corresponding  confi- 


Fig.  1.  Condition  monitoring  structure  using  Hotelling  method. 
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dence  level)  for  conformity  or  to  eliminate  the  predictions  with 
T 2  statistics  exceeding  the  upper  control  limit.  At  phase  two, 
when  an  online  measurement  of  the  output  voltage  is  obtained 
for  monitoring,  it  will  be  incorporated  into  the  Hotelling  T 2  anal¬ 
ysis  to  establish  a  revised  control  limit  that  takes  into  account 
the  increase  of  dimension  of  the  dataset  due  to  the  inclusion 
of  the  new  measurement.  Statistically,  the  Hotelling  T 2  analysis 
enables  the  conformity  check  between  the  voltage  measurement 
and  the  baseline  healthy  voltage  prediction.  Fault  condition  will 
then  be  declared  if  the  T 2  statistics  of  the  online  voltage  measure¬ 
ment  of  the  fuel  cell  being  monitored  exceeds  the  upper  control 
limit.  Clearly,  this  proposed  condition  monitoring  system  can 
deal  with  the  two  major  issues  in  fuel  cell  applications,  the  com¬ 
plex  multi-physics  operation  and  the  system  uncertainties.  The 
components  shown  in  Fig.  1  will  be  explained  in  detail  in  what 
follows. 

2.2.  Reference  model  development 

Without  loss  of  generality,  we  use  the  Ballard  Mark  V  5  kW 
stack  as  benchmark  example  [13]  and  develop  the  mathematical 
model.  The  reason  we  choose  this  particular  stack  for  demon¬ 
strating  the  methodology  development  is  that  the  relevant  param¬ 
eters  and  typical  operating  processes  are  well-documented.  The 
stack  is  composed  of  35  cells.  In  a  previous  study,  a  general 
system-level  modeling  methodology  for  PEM  fuel  cell  has  been 
outlined  [14].  In  this  section,  we  develop  a  mathematical  model 
that  is  specifically  tailored  for  stack  health  monitoring.  It  has 
been  pointed  out  that,  when  only  the  stack  voltage  is  used 
as  information  carrier,  it  is  difficult  to  detect  the  failure  of  a 
single  cell  in  a  large  stack  [4].  The  reason  is  that  the  wiring 
could  become  extremely  complicated  if  one  wants  to  measure 
single-cell  output  voltage.  It  is  suggested  that  the  group  volt¬ 
age  monitoring  of  no  more  than  five  cells  seems  to  be  feasible 
from  stack  operation  standpoint  [4].  Therefore,  in  our  study  the 
stack  modeling  will  be  based  on  the  group  consideration  of  five 
cells.  The  dry  fuel  from  a  compressed  gas  tank  and  the  dry  puri¬ 
fied  air  from  an  air  compressor  enter  the  stack  separately  and 
flow  directly  into  the  humidification  section  where  the  fuel/gas 
streams  pick  up  water.  The  fuel/gas  main  streams  are  further 
separated  and  then  individually  fed  into  each  of  the  fuel  cells  in 
their  respective  anode  and  cathode  channels  in  the  active  fuel  cell 
section.  The  electrochemical  reaction  takes  place  at  the  individ¬ 
ual  MEA.  The  electrical  power  is  then  generated,  and  the  surplus 
fuel/gas  and  produced  water  due  to  the  electrochemical  reaction 
flow  out  through  their  respective  outlets. 

PEM  fuel  cell  has  been  modeled  at  various  levels  of  complex¬ 
ity  [1-3,14-17].  While  the  micro-level  models  of  single  cell  aim 
at  the  cell  component  synthesis/design,  these  models  consist  of 
complex  partial  differential  equations  and  are  mostly  based  on 
steady  state  conditions.  Moreover,  the  model  complexity  dra¬ 
matically  increases  after  one  packs  up  the  single-cell  models 
into  the  stack  model.  These  models  usually  require  very  high 
computational  cost  and  cannot  be  used  in  online  monitoring  that 
needs  real-time  prediction.  In  order  to  establish  a  tractable  stack 
model  that  can  be  directly  used  in  health  monitoring,  here  we  use 
the  control  volume  (CV)  approach  at  the  modeling  stage.  The 


aforementioned  benchmark  stack  is  divided  into  seven  groups, 
each  of  which  is  composed  of  five  cells.  The  fuel  cell  sub-group  is 
divided  into  three  control  volumes,  the  anode  channel,  the  cath¬ 
ode  channel,  and  the  fuel  cell  body.  As  each  fuel  cell  sub-group 
is  composed  of  five  single  fuel  cells,  the  five  anode  channels 
involved  in  this  sub-group  are  treated  as  one  anode  channel  CV. 
Same  is  applied  to  the  cathode  channel  CV  and  the  fuel  cell 
body  CV.  For  a  sub-group  composed  of  five  single  fuel  cells, 
the  species  conservation  and  energy  conservation  principles  are 
applied  to  three  CVs  [14].  If  pure  hydrogen  is  used  as  fuel  and  is 
humidified  before  entering  the  anode  channel,  two  species,  i.e., 
hydrogen  and  water  (liquid  and  vapor),  will  be  involved  in  the 
modeling.  There  are  three  equations  for  mass  conservation  and 
one  equation  for  energy  conservation.  It  is  worth  mentioning  that 
phase  change  plays  a  critical  role  in  PEM  fuel  cell  operations. 
This  is  especially  relevant  in  a  model-based  condition  modeling 
because  the  water  content  and  heat  directly  decide  the  proper 
operation  of  a  PEM  fuel  cell.  The  present  model  explicitly  takes 
this  into  account. 


2.2.7.  Anode  channel 
Mass  conservation: 


dmn2 
d  t 


NI 

Wn2,out  -  — 


dm 


an,v 


dt 


—  Wan, v, in  Wan^ont  Wv?mem  T  {sgn}p  Wan  V  phase 


dmanw 

dt 


{s£nlp  Wan,  v,  phase  Wan,w  ,out 


All  notations  are  defined  in  the  Nomenclature.  In  the  above 
equations,  the  hydrogen  inlet  flow  rate  WH2,m  depends  on  the 
source  fuel  and  anode  channel  pressure,  relative  humidity,  tem¬ 
perature,  and  anode  inlet  flow  rate  coefficient.  The  flow  rates 
of  other  species,  WH2,0ut,  Wan,V5in,  Wan,v,out,  and  Wan,w,out,  can 
be  obtained  similar  to  WH2,in-  Wv,mem  is  the  water  flow  rate 
across  the  membrane  that  can  be  calculated  based  on  the  electro- 
osmotic  drag  force  and  the  back-diffusion  from  the  cathode  side 
to  the  anode  side  caused  by  the  water  content  difference  [16,17]. 

In  Eqs.  (2)  and  (3),  the  sign  function  {sgn}p  is  used  to  charac¬ 
terize  the  direction  of  phase  change,  (+1)  indicating  evaporation 
and  (—1)  indicating  condensation.  It  is  worth  mentioning  that 
the  two-phase  change  is  an  important  feature  in  PEM  fuel  cells 
that  operate  at  relatively  low  temperatures.  It  plays  an  impor¬ 
tant  role  in  heat  and  water  management,  and  humidification  of 
membrane  and  electro-catalyst  layers  [14].  In  literature,  various 
modeling  approaches  have  been  proposed  to  characterize  this 
phenomenon.  In  microscopic  level  modeling,  multi-phase  mix¬ 
ture  based  two-phase  flow  model  is  employed,  where  the  phases 
are  assumed  to  be  distinct  and  separable  components  with  non¬ 
zero  interfacial  areas,  and  the  phase  change  is  characterized  by 
smoothly  varying  phase  compositions  in  their  mixture  [18,19]. 
Such  micro-level  models  generally  reflect  the  local  behavior 
only  and  require  very  high  computational  cost  which  is  not 
feasible  for  online  applications.  There  has  been  recent  interest 
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in  developing  system-level  models  that  can  provide  the  global 
characterization  of  the  system  behavior  [14,16].  Pukrushpan  et 
al.  [16]  developed  a  constant  temperature  model  where  it  was 
assumed  that  the  vapor  condenses  into  the  liquid  form  once  the 
relative  humidity  of  the  gas  exceeds  100%.  Xue  and  Tang  [14] 
proposed  a  phase  change  modeling  approach,  which  is  incorpo¬ 
rated  into  a  comprehensive  system-level  dynamic  model  of  PEM 
fuel  cells.  The  fundamental  principle  is  that  if  the  vapor  partial 
pressure  is  larger  than  the  in  situ  saturation  vapor  pressure,  liquid 
water  will  be  formed  through  condensation.  On  the  other  hand, 
if  the  vapor  partial  pressure  is  smaller  than  the  corresponding  in 
situ  saturation  pressure,  the  existent  liquid  water  will  evaporate 
until  the  vapor  partial  pressure  reaches  the  saturation  pressure. 
In  this  study,  we  employ  this  modeling  approach,  and  the  phase 
change  rate  equation  is  given  as  [14] 


Wi 


an,  v,  phase 


tfv Tan  j(Ti 


an? 


1 

7an,v> 


_ /  pTm 

\  s1  an,  v, sat 


where  Van  is  the  channel  volume  of  the  sub-group. 
Energy  conservation: 
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+  (AA)an(7’body-  7^)-  Wv  ,mem%),v  [{sgn}dran 
+  ({sgn}d  —  l)rt,0dy] 


where  the  sign  function  {sgnjd  is  used  to  indicate  the  vapor 
transport  direction  in  MEA,  {sgnjd  =  +1  if  vapor  is  transported 
from  the  anode  channel  to  the  cathode  channel,  and  {sgnjd  =  0 
otherwise.  (M)ca  represents  the  convective  heat  transfer  coeffi¬ 
cient  between  the  cathode  channel  control  volume  and  the  fuel 
cell  body. 


2.2.2.  Cathode  channel 
Mass  conservation: 
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Energy  conservation 
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2.2.3.  Fuel  cell  body 

The  fuel  cell  body  is  treated  as  a  solid  control  volume.  The 
mass  of  the  fuel  cell  body  is  assumed  to  be  a  constant,  and  thus 
the  mass  conservation  is  trivial.  The  only  concerned  dynamic 
behavior  of  fuel  cell  body  is  the  temperature.  Using  the  energy 
conservation,  we  have 


d7b0dy  NI  NI 
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NI 
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T  Tfjody )  +  (hA)cw(Tcw  Tdody) 

+  (KA)(r'-d1y  -  2r'ody  +  r'+‘y)  (9) 


If  the  water  is  transported  from  the  anode  side  to  the  cath¬ 
ode  side,  i.e.,  {sgnjan=  1,  the  anode  channel  temperature  has 
contribution  to  the  fuel  cell  temperature  change.  Otherwise, 
{sgn}an  =  0,  and  it  is  the  cathode  channel  temperature  rca  that 
will  affect  the  fuel  cell  temperature.  The  last  term  in  Eq.  (9) 
represents  the  conduction  heat  transfer  among  the  neighboring 
sub-groups. 

2.2.4.  Electrochemical  reaction 

Typically,  the  fuel  cell  output  voltage  is  the  summation  of 
three  effects  [16],  the  Nernst  potential,  the  cathode  and  anode 
activation  overvoltage,  and  the  ohmic  overvoltage  due  to  internal 
resistance, 

Ecell  —  ^Nernst  —  ?7act  —  7?ohmic  (10) 

where  the  Nernst  potential,  the  activation  voltage  and  the  ohmic 
voltage  loss  are  functions  of  the  operating  conditions.  One  way 
of  modeling  the  activation  overvoltage  and  ohmic  overvoltage  is 
to  take  into  account  that  the  membrane  conductivity  is  strongly 
dependent  upon  the  membrane  water  content  and  fuel  cell  tem¬ 
perature  [20].  Alternatively,  the  expressions  of  the  activation 
overvoltage  and  the  ohmic  overvoltage  may  be  characterized 
using  the  experimental  data  under  various  cell  temperature  and 
load  current  [13,14,16], 

7? act  —  —  §1  —  §27body  —  §3  Tbody  ln(c02 )  —  §4  Tbody  (7)  (H) 

7? ohmic  —  ms  3“  ?6  Tbody  T  §7  7)  (12) 

where  Cq2  is  the  concentration  of  oxygen  at  the  catalyst  interface, 
and  §/  (/=  1,  . . .,  7)  are  parameters  calculated  from  the  experi¬ 
mental  data  using  linear  regression.  It  should  be  noted  that  the 
electrochemical  effects  may  vary  slightly  from  cell  to  cell  and 
depend  on  the  complicated  operating  conditions  within  a  cell. 
In  order  to  take  into  account  the  system  variation/uncertainty,  in 
this  research  we  use  the  linear  regression  expressions  where  the 
coefficients  can  be  easily  perturbed  during  the  statistical  sam- 
pling. 
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As  the  cells  are  serially  connected,  the  sub-group  output  volt¬ 
age  will  be  the  summation  of  the  five  single  cell  output  voltages. 

2.2.5.  Stack  model  assembly 

Eqs.  (1)-(12)  provide  a  basis  for  the  modeling  of  the  con¬ 
cerned  fuel  cell  sub-group.  The  assembly  of  such  sub-group 
systems  will  form  the  stack  model,  which  also  include  heat  con¬ 
duction  effect  among  these  fuel  cell  sub-groups.  As  the  total 
inlet  flow  rate  of  fuel/gas  is  individually  fed  into  each  of  the 
fuel  cells  and  the  exhaustion  of  each  fuel  cell  converges  to  the 
stack  outlet,  the  inlet/outlet  flow  rates  of  the  stack  will  be  the 
summation  of  the  respective  sub-group  inlet/outlet  flow  rates. 
The  cooling  water  flow  rate  is  treated  in  the  same  manner.  We 
then  have, 
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In  all,  there  are  seventy  first-order  dynamic  equations  for  the 
stack.  A  variety  of  stack  state  variables  can  be  obtained  from  the 
model  simulation.  In  this  study,  the  sub-group  output  voltage 
will  be  employed  to  monitor  the  conditions  of  the  stack  using 
the  Hotelling  statistical  analysis. 

2.3.  Baseline  simulation  with  system  uncertainties 

In  the  aforementioned  dynamic  model,  all  the  system  param¬ 
eters  are  assumed  as  deterministic  parameters.  However,  in  a  real 
fuel  cell  system,  variations/uncertainties  exist.  Part  of  the  uncer¬ 
tainties  comes  from  the  normal  time-evolution  nature  of  fuel  cell 
components  as  well  as  the  cell-to-cell  manufacturing  variation. 
From  the  modeling  perspective,  even  more  significant  uncertain¬ 
ties  may  appear  since  the  fuel  cell  model  involves  coefficients 
obtained  based  on  experimental  studies  (e.g.,  in  Section  2.2.4). 
PEM  fuel  cells  are  highly  non-linear.  The  practical  operation 
of  a  PEM  fuel  cell  is  essentially  dynamic  and  thus  the  operat¬ 
ing  conditions  cover  various  combinations  of  parameters/states 
that  may  deviate  from  the  nominal  condition.  In  addition,  all  the 
sensor  measurements  (of  the  input  state  variables  and  the  output 
voltage)  may  be  subjected  to  noise. 

In  order  to  take  into  account  the  fuel  cell  varia¬ 
tions/uncertainties  effect  in  the  condition  monitoring  system, 
in  this  research  we  let  the  key  empirical  parameters  used  in 
the  baseline  model  be  random  variables  with  specified  statistics 


(Table  2).  As  a  result,  the  model  outputs  become  stochastic  in 
nature.  Stochastic  analysis  is  a  procedure  where  the  input  para¬ 
metric  uncertainties  propagate  through  a  physical  process  (such 
as  the  fuel  cell  system).  As  most  of  the  physical  processes  in 
the  PEM  fuel  cell  system  are  governed  by  differential/partial 
differential  equations  and  are  non-linear  processes,  the  uncer¬ 
tainty  propagation  using  the  closed- form  result  is  not  feasible. 
In  the  proposed  condition  monitoring  system,  we  will  use  the 
sampling  approach  to  perturb  the  uncertain  parameters,  which 
essentially  leads  to  a  group  of  output  voltages  as  the  baseline. 
This  group  of  output  voltages  covers  the  possible  output  range 
under  the  normal  operation  condition,  which  will  be  used  in  the 
Hotelling  T 2  statistics  analysis  for  damage  detection. 

In  traditional  stochastic  analysis,  the  samples  are  selected 
randomly  using  such  as  the  Monte  Carlo  technique  [21].  The 
computational  cost  for  such  sampling,  however,  is  in  general 
quite  high  [21-23] .  To  improve  the  sampling  efficiency,  the  strat¬ 
ified  sampling  methods  have  been  developed,  such  as  the  Latin 
hypercube  sampling  (LHS)  method  [24]  and  the  Hammersley 
sequence  sampling  (HSS)  method  [25].  The  LHS  method  is  a 
stratified  sampling  method,  where  the  distribution  is  divided  into 
m  intervals  with  equal  probability,  and  the  sample  is  then  picked 
randomly  from  each  interval.  In  the  HSS  method,  m  sampling 
points  are  uniformly  placed  in  a  k-dimensional  cube  through  a 
low-discrepancy  design,  and  the  representative  sample  is  gener¬ 
ated.  In  order  to  demonstrate  the  proposed  condition  monitoring 
strategy,  here  we  assume  that  the  empirical  parameters  listed  in 
Table  2  all  have  the  Gaussian  distribution,  and  the  corresponding 
standard  deviations  are  also  listed  in  the  same  table.  Let  the  total 
number  of  uncertain  parameters  be  q.  Using  the  LHS  procedure, 
m  values  of  each  of  the  parameters  are  sampled  under  the  Gaus¬ 
sian  distribution  assumption.  We  then  obtain  a  m  x  q  sampling 
data  matrix.  For  each  column  data,  the  m  values  are  randomly 
distributed  with  one  from  each  interval  as  mentioned  above,  and 
they  are  randomly  permuted.  When  a  set  of  the  input  measure¬ 
ment  data  (from  physical  fuel  cell  system)  comes  in,  the  system 
model  will  be  simulated  m  times  using  each  row  of  the  above 
m  x  q  sampling  data  matrix.  Correspondingly,  we  will  obtain  m 
output  voltages,  which  will  be  used  as  the  baseline  for  damage 
detection.  Clearly,  these  m  output  data  will  cover  the  uncertainty 
band  resulted  from  the  system  parameter  uncertainty  propaga¬ 
tion  through  the  system  model,  and  the  computational  time  is 
reduced  using  the  LHS  procedure  for  sampling. 

2.4.  Hotelling  T2  statistical  analysis  for  stack  condition 
monitoring 

The  multiple  simulations  of  PEM  stack  model  under  statis¬ 
tical  sampling  of  parametric  uncertainties  will  yield  multiple 
output  voltage  time-responses  as  the  baseline.  These  voltage 
responses  (as  functions  of  time,  due  to  the  dynamic  nature  of 
the  PEM  fuel  cell)  will  then  be  compared  with  the  real-time 
measurement  of  the  output  voltage  for  condition  monitoring. 
Hotelling  first  proposed  the  multi-variate  statistical  analysis  for 
conformity  check  [10].  This  method  is  considered  robust  and 
highly  sensitive  for  such  task,  and  has  been  widely  used  in  pro¬ 
cess  control  and  fault  detection  [11,12].  The  analysis  is  based 
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upon  the  T 2  statistics  and  the  comparison  with  certain  control 
limits. 

In  general,  for  a p-  variable  problem  with  each  variable  having 
m  observations,  the  mean  vector  is  defined  as, 

x  =  [xi,x2,  ...,x/)]T  (16) 

where  xj  =  (1  / m)Yl?=ixij  is  the  mean  estimation  for  the  j- th 
variable.  The  sample  covariance  matrix  S  can  be  calculated  as, 

I  m 

S  = - Ty'(x,-  -  x)(Xj  -  x)T  (17) 

m  —  1 

i=l 

The  covariance  matrix  S  indicates  the  relationship  existing 
among  the  p- variables.  In  the  present  study,  the  output  voltage 
of  each  sub-group  is  used  as  the  information  carrier  for  condi¬ 
tion  monitoring.  Therefore,  the  output  voltage  history  calculated 
or  measured/sampled  will  be  grouped,  with  each  group  having 
p  time  points.  We  will  perform  condition  monitoring  based  on 
these  grouped  p- variables  successively.  Recall  that  in  order  to 
take  into  account  the  system  uncertainties  caused  by  the  empir¬ 
ical  modeling,  the  fuel  cell  model  will  be  simulated  multiple 
times  with  statistical  sampling  (through  the  LHS  procedure)  of 
physical  parameters.  We  may  then  obtain  m  voltage  output  obser¬ 
vations  correspondingly.  There  are  two  distinct  phases  involved 
in  constructing  the  Hotelling  T 2  control  limits  [10].  First,  we  use 
the  aforementioned  m  observations  (simulations)  of  healthy  fuel 
cell  model  to  establish  a  Phase  I  baseline  and  an  upper  control 
limit  under  a  specified  confidence  level.  The  purpose  is  to  exam¬ 
ine  the  baseline  data  and  also  the  relation  of  normal  condition 
with  respect  to  the  upper  control  limit.  For  a  given  voltage  output 
simulation/observation  x/,  the  statistic  T 2  is  defined  as 

T2  =  (xi  -  x)TS-\Xi  -  x)  (18) 

It  can  be  proved  that  the  statistic  T 2  follows  the  F-distribution. 
Therefore,  under  a  given  Type  1  error  probability  a ,  the  Phase  I 
upper  control  limit  UCLlp  m  a  can  be  established  as  [10] 


Table  1 

Parameters  used  in  the  model  validation  and  simulation 


Parameters 

Value 

Reference 

Fin, in  (kg(Pas)-1) 

1.23  X  10“10 

a 

Psan  (psig) 

35 

[13] 

Fn, in  (°C) 

23.5 

[13] 

ty^an 

0.9 

b 

Fa, out  (kg(Pas)-1) 

2  x  nr5 

a 

Tc a, in  (°C) 

23.5 

[13] 

<Pca 

0.0 

b 

Wcw  (kgs-1) 

0.0331568 

[13] 

Wcw  (kg) 

2.5 

b 

(MJanCWK-1) 

2 

[13] 

(M)oo(WK-1) 

17 

[13] 

Too  (°C) 

23.5 

[13] 

t,n  (m) 

1.275  x  nr4 

[16] 

Pm, dry  (kgcm-3) 

0.002 

[1] 

Fn,out  (kg(Pas)-1) 

2  x  nr5 

a 

Tatm  (Pd) 

1.013  x  105 

[13] 

Fn  (m3) 

0.005 

[16] 

Fa, in  (kg(Pas)-1) 

9.148  x  10“10 

a 

Psca  (psig) 

35 

[13] 

Vca  (m3) 

0.01 

[16] 

MC  (JK-1) 

35000 

[13] 

Tew, in  (°C) 

23.5 

[13] 

N 

35 

[13] 

(M)ca( WK->) 

10 

[13] 

(M)cw(WK->) 

50 

[13] 

T  (S) 

30 

b 

Fell  (m2) 

0.0232 

[13] 

Mm, dry  (kgmol-1) 

1.1 

[1] 

a  Coefficients  are  so  chosen  that  the  fuel/gas  inlet  flow  rate  is  the  same  as  those 
of  the  experimental  setup. 
b  Parameters  are  assumed. 


the  analysis.  Since  the  new  voltage  measurement  is  indepen¬ 
dent  of  the  Phase  I  baseline,  the  Phase  II  upper  control  limit  is 
modified  as  [10] 


UCL 


2 

p,m,a 


p(m  -  1  )(m  +  1) 

- - - - - Fa(p,  m-  p) 

m(m  —  p) 


(20) 


i  p(m  —  1) 

UC \2  =  Fa(p ,  m-p )  (19) 

1  ’  m(m  —  p) 

where  Fa(p ,  m—p)  is  the  1  —  a  percentile  of  the  F-distribution 
with  p  and  m—p  degrees  of  freedom.  Indeed,  we  may  use  m 
observations  of  the  output  voltage  under  normal  operating  con¬ 
dition  predicted  by  the  fuel  cell  model  to  establish  a  control  limit 
with  1  —  a  confidence  level.  At  this  Phase  I  stage,  the  observa¬ 
tions  of  normal  fuel  cell  output  (obtained  by  model  prediction 
under  statistical  sampling)  will  undergo  a  self-checking  proce¬ 
dure.  The  statistic  T 2  of  every  output  voltage  prediction  will  be 
compared  with  the  upper  control  limit.  If  the  statistic  T2  of  cer¬ 
tain  prediction  exceeds  the  upper  control  limit,  we  may  either 
eliminate  this  prediction  from  the  baseline  data  or  revise  the 
upper  control  limit  (and,  correspondingly,  the  confidence  level). 
This  procedure  iterates  until  the  baseline  data  are  purified. 

The  baseline  output  voltages  under  normal  operating  condi¬ 
tion  will  then  be  utilized  to  establish  Phase  II  statistic  T 2  for 
condition  monitoring.  This  is  facilitated  by  incorporating  the 
real-time  online  measurement  of  fuel  cell  output  voltage  into 


In  other  words,  when  a  real-time  measurement  of  output  volt¬ 
age  comes  in,  we  need  to  use  Eq.  (18)  to  calculate  the  statistic 
T2  value  of  this  new  p-  variable.  If  this  value  exceeds  the  Phase  II 
upper  control  limit  that  we  have  established  in  Eq.  (20),  we  may 
conclude,  with  1  —  a  confidence  level,  that  the  voltage  output  is 
abnormal  and  fault  condition  has  occurred  in  the  corresponding 
sub-group  of  the  stack.  From  the  above  procedure,  we  can  see 
that  the  Hotelling  T2  analysis  allows  the  direct  comparison  of 


Table  2 

Parameters  used  in  the  simulation  and  uncertainties 


Parameters 

Value 

Standard 
deviation  (%) 

Reference 

0.944 

0.01 

[13]  for  mean  value 

£2 

-3.54  X  10“3 

0.1 

[13]  for  mean  value 

£3 

-4.68  x  10“4 

0.1 

[13]  for  mean  value 

£4 

1.96  x  10“4 

0.1 

[13]  for  mean  value 

£5 

-3.3  x  10“3 

0.1 

[13]  for  mean  value 

^6 

7.55  x  10"6 

0.1 

[13]  for  mean  value 

£7 

-1.1  x  10"6 

0.1 

[13]  for  mean  value 
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the  online  measurement  with  a  group  of  possible  output  voltage 
under  normal  operation  conditions,  which  is  a  key  feature  of  the 
proposed  condition  monitoring  strategy. 

3.  Simulation  and  discussion 

3.1.  Dynamic  simulation  of  fuel  cell  stack  using  nominal 
parameters 

In  order  to  explore  the  dynamic  behavior  of  the  PEM  fuel  cell 
stack  under  various  fault  conditions,  a  series  of  step  changes  in 
load  current  is  applied  to  the  fuel  cell  system  to  emulate  the  exter¬ 
nal  load  change.  The  initial  conditions  and  system  parameters 
are  listed  in  Tables  1  and  2.  Here  in  this  first  study,  we  assume  an 
ideal  system  without  uncertainties,  and  hence  all  system  param¬ 
eters  take  the  nominal  values.  The  purpose  is  to  compare  the 
nominal  difference  of  fuel  cell  state  variables  caused  by  the  fault 
conditions.  During  the  simulation,  the  fault  condition  is  delib¬ 
erately  introduced  into  the  system.  The  total  time  period  for 
the  simulation  is  30,000  s.  Initially,  no  external  load  current  is 
applied  to  the  system.  At  the  100th  s,  the  system  starts  up  with  a 
20  A  external  load  current.  At  the  5000th  s,  the  membrane  drying 
of  a  single  cell  is  introduced  into  sub-groups  3  and  6,  respec¬ 
tively.  This  fault  condition  recovers  to  the  normal  condition  at 
the  7000th  s.  At  the  1 1,000th  s,  the  load  current  continuously 
rises  up  to  30  A.  The  cathode  flooding  is  then  introduced  to  a 
single  cell  in  sub-group  6  at  the  17,000th  s  which  continues  for 
3000  s.  The  fuel  cell  then  returns  to  the  normal  condition.  The 
last  fault  introduced  is  the  cathode  flooding  of  a  single  cell  of 
sub-group  3  at  the  25,000th  s.  A  number  of  stack  state  variable 
histories  are  shown  in  Figs.  2-13. 

We  first  analyze  the  fuel  cell  state  variable  histories  to  gain  the 
physical  insight.  Fig.  2  shows  the  stack  temperature  history.  All 
seven  sub-groups  show  similar  trend,  i.e.,  the  stack  temperature 
increases  following  the  increase  of  the  external  load  current.  The 
temperature  is  influenced  by  the  occurrence  of  cell  faults.  At  the 
5000th  s,  the  membrane  drying  is  introduced  into  sub-groups  3 
and  6,  respectively,  and  the  drying  level  of  sub-group  3  is  higher 
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310 


0 

s_ 

0 

0 

l_ 

0 

Q. 

E 

0 


Time(sec) 


-©- 

Group  1 

— 

Group  2 

Group  3 
Group  4 

-ft- 

Group  5 

-A- 

Group  6 

-v- 

Group  7 

J _ I 

2  2.5 


3 

x  104 


x  104 


1.5 

Time(sec) 


2.5  3 

x  104 


Fig.  2.  Stack  temperature  history  simulation. 


Fig.  5.  Oxygen  partial  pressure  simulation  in  cathode  CV. 
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Fig.  6.  Vapor  partial  pressure  simulation  in  cathode  CV. 


Fig.  9.  Direct  comparison  of  sub-group  voltage  values. 
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Fig.  10.  Baseline  self-checking  result  of  sub-group  3  voltage. 
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Fig.  11.  Baseline  self-checking  result  of  sub-group  6  voltage. 
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Fig.  12.  Hotelling  analysis  of  sub-group  3  voltage. 


of  the  consumed  fuel  and  water  generated  due  to  the  electro¬ 
chemical  reaction.  Consequently,  the  hydrogen  partial  pressure 
in  the  anode  CV  decreases  even  though  the  anode  CV  tempera¬ 
ture  increases,  while  the  partial  pressure  of  the  vapor  increases 
due  to  both  anode  CV  temperature  increase  and  the  increase  of 
water  diffusion  effect  from  the  cathode  side  to  the  anode  side. 
In  Figs.  5  and  6,  the  oxygen  and  vapor  partial  pressures  in  the 
cathode  show  similar  trends  to  those  of  the  anode  CV. 

Figs.  7  and  8  illustrate  the  liquid  water  generated  due  to  phase 
change  in  the  anode  and  cathode  CVs,  respectively.  It  is  apparent 
that  there  is  no  liquid  water  in  the  anode  CV,  while  the  liq¬ 
uid  water  generated  in  the  cathode  CV  increases  following  the 
increase  of  the  external  load  current.  Clearly,  as  the  electrochem¬ 
ical  reaction  intensifies,  the  water  vapor  generated  increases  and 
then  the  condensation  takes  place  in  the  cathode  CV.  One  may 
see  that  the  faults  are  explicitly  simulated  using  the  above  devel¬ 
oped  model,  and  they  do  cause  some  interesting  changes  during 
fuel  cell  operations. 


than  that  of  the  sub-group  6.  It  is  apparent  that  temperatures 
of  sub-groups  3  and  6  both  increase  with  sub-group  3  having 
a  higher  temperature.  Similarly,  the  flooding  of  cathode  elec¬ 
trode  blocks  the  transport  of  oxygen  from  the  cathode  channel 
to  the  cathode  electrolyte  and,  consequently,  the  oxygen  partial 
pressure  at  the  electrolyte  decreases,  which  directly  leads  to  the 
output  voltage  decrease.  Here  it  is  assumed  that  the  oxygen  at  the 
electrolyte  is  adequate  for  the  electrochemical  reaction  under  the 
loading  conditions.  The  decrease  of  the  electrical  power  leads 
to  the  temperature  increase  of  sub-groups  3  and  6,  provided  that 
the  consumed  fuel  is  the  same  as  that  of  the  normal  operating 
conditions.  This  result  is  shown  in  Fig.  2  at  the  17,000th  and 
25,000th  s,  respectively.  One  may  observe  that  the  temperature 
increase  of  sub-groups  3  and  6  influences  the  neighboring  stack 
temperature. 

Figs.  3  and  4  show  the  hydrogen  and  vapor  partial  pressures  in 
the  anode  channel  CV  of  seven  sub-groups.  The  hydrogen  partial 
pressure  tends  to  go  down  following  the  load  current  increase. 
The  increase  of  the  load  current  directly  leads  to  the  increase 


Fig.  13.  Hotelling  analysis  of  sub-group  6  voltage. 


3.2.  Hotelling  monitoring 

In  this  sub-section,  we  demonstrate  the  Hotelling  statistical 
analysis  method  for  fault  detection.  As  stated,  under  the  respec¬ 
tive  fault  conditions  in  three  time  periods  (5000th-7000th  s,  cells 
in  sub-groups  3  and  6  having  faults;  17,000th-20,000th  s,  a  cell 
in  sub-group  6  having  fault;  25, 000th-28, 000th  s,  a  cell  in  sub¬ 
group  3  having  fault),  the  output  voltages  of  sub-groups  have 
very  small  fault-caused  variations  relative  to  the  healthy  con¬ 
dition  (the  exact  maximum  variation  using  nominal  parameters 
without  uncertainty  is  0.0446  V).  In  practical  applications,  the 
healthy  responses  are  subjected  to  model  uncertainty,  and  both 
the  healthy  and  faulty  responses  will  be  contaminated  by  mea¬ 
surement  noise.  In  this  research,  the  modeling  uncertainty  is 
accounted  for  by  assuming  variance  in  the  empirical  parameters 
used  in  the  description  of  electrochemical  reaction  (Table  2).  It 
should  be  noted  that  little  knowledge  currently  exists  in  quanti¬ 
fying  such  uncertainty,  and  the  standard  deviation  values  shown 
in  Table  2  are  obtained  by  trial  and  error  that  yield  reason¬ 
able  output  voltage  history  calculation  in  the  simulation.  In 
future,  comprehensive  experiments  are  needed  in  quantifying 
such  uncertainty.  In  addition  to  such  parametric  uncertainties, 
in  this  study  we  assume  that  the  voltage  measurement  is  sub¬ 
jected  to  random  noise  with  0.016  V  standard  deviation.  Fig.  9 
illustrates  the  output  voltages  of  the  faulty  sub-groups  3  and  6, 
respectively.  Also  shown  are  the  healthy  output  voltages  selected 
from  the  corresponding  baseline  (obtained  via  statistical  sam¬ 
pling,  which  will  be  further  discussed  later)  under  the  same 
operating  conditions.  Clearly,  the  fault  conditions  cause  very 
insignificant  change  of  output  voltage,  which  is  further  con¬ 
taminated  by  the  noise.  Therefore,  it  is  virtually  impossible  to 
declare  fault  occurrence  based  upon  the  direct  comparison  of 
output  voltages,  as  shown  in  Fig.  9. 

We  then  implement  the  Hotelling  statistic  T2  detection 
method,  where  we  use  sampled  output  voltage  measurements  as 
the  p-  variable.  Here,  for  easy  identification  of  the  fault  sources, 
the  Hotelling  detection  is  applied  for  each  sub-group  of  the  volt¬ 
ages.  In  order  to  establish  the  baseline  and  the  control  limit, 


398 


X.  Xue  et  al.  /  Journal  of  Power  Sources  162  (2006)  388-399 


a  set  of  fuel  cell  output  voltages  are  generated  using  the  fuel 
cell  model  under  statistical  sampling  of  parametric  uncertain¬ 
ties  (through  the  LHS  procedure)  under  the  normal  operating 
conditions.  Here,  the  Gaussian  distributions  are  assumed  for  the 
uncertain  parameters  (shown  in  Table  2)  and  each  distribution  is 
divided  into  16  strata.  Output  voltages  are  grouped  every  6  time 
points,  i.e.,  p-  6.  Hotelling  analysis  are  performed  successively 
for  the  entire  voltage  histories.  After  self-checking,  a  baseline 
composed  of  16  sample  voltages  is  constructed.  Figs.  10  and  1 1 
show  the  corresponding  self-checking  results  of  normal  voltage 
outputs  of  sub-groups  3  and  6,  respectively.  Here,  the  confi¬ 
dential  level  of  the  F-distribution  is  set  to  be  95%.  The  degree 
of  freedom  of  the  Phase  I  self-checking  is  p  x  m  =  6  x  16.  The 
upper  control  limit  of  Phase  I  is  27. 14.  Clearly,  the  self-checking 
results  of  the  baseline  are  below  the  Phase  I  upper  control  limit, 
and  we  may  now  declare  with  95%  confidence  that  the  baseline 
is  purified. 

At  Phase  II,  the  (simulated)  online  measurement  of  sub-group 
output  voltage  measurements  is  inserted  into  the  Hotelling  anal¬ 
ysis  module,  and  analyzed  based  on  the  baseline  developed  at 
Phase-I  stage.  If  the  statistic  T2  of  the  online  output  voltage 
measurement  exceeds  the  Phase  II  upper  control  limit,  a  fault 
condition  can  be  declared;  otherwise,  the  operating  condition 
is  treated  as  normal.  Again,  each  subgroup  output  voltage  of 
the  online  fuel  cell  output  voltage  is  consecutively  grouped  for 
every  6  time  points.  One  advantage  of  this  strategy  is  that  the 
condition  monitoring  system  would  have  enough  time  to  pro¬ 
cess  all  the  calculations.  With  the  new  online  measurement,  the 
Phase  II  upper  control  limit  can  be  calculated  as  30.76  under 
a  95%  confidence  level  (Eq.  (20)).  We  then  test  this  condition 
monitoring  method  by  incorporating  the  simulated  voltage  mea¬ 
surement  of  the  faulty  fuel  cell  system.  Figs.  12  and  13  show  the 
Hotelling  results  of  abnormal  voltage  outputs  of  sub-groups  3 
and  6,  respectively.  It  can  be  clearly  seen  that  the  T 2  values  of  the 
measurements  exceed  the  upper  control  limit  at  the  time  instant 
when  the  fault  occurs.  The  Hotelling  statistic  T 2  of  sub-groups  3 
and  6  declare  both  fault  at  the  5000th  s.  The  fault  condition  lasts 
for  about  2000  s  and  then  the  T 2  values  drop  to  below  the  upper 
control  limit.  At  the  17,000th  s,  the  Hotelling  statistic  T2  of  sub¬ 
group  6  goes  beyond  the  upper  control  limit,  which  lasts  for 
3000  s.  At  the  25,000th  s,  the  Hotelling  statistic  T2  of  sub-group 
3  declares  fault  condition.  While  the  direct  comparison  of  the 
fault  and  normal  voltages  in  Fig.  9  cannot  give  definitive  conclu¬ 
sion,  the  Hotelling  method  significantly  improves  the  detection 
sensitivity  and  can  correctly  predict  the  fault  occurrence  imme¬ 
diately,  even  under  system  parameter  uncertainties  and  noisy 
measurements. 

4.  Concluding  remarks 

This  research  proposes  a  condition  monitoring  approach  for 
the  PEM  fuel  cells.  The  basic  idea  is  to  compare  the  fuel  cell  out¬ 
put  voltage  measurement  with  model  prediction  under  the  same 
operating  conditions.  We  first  develop  a  fuel  cell  stack  model 
which  can  simulate  the  complicated  transient  behavior  and 
dynamic  interaction  of  the  temperature,  gas  flow,  phase  change 
in  the  anode  and  cathode  channels,  and  membrane  humidifica¬ 


tion  under  operating  conditions.  Using  this  model  as  basis,  we 
then  simulate  the  stack  output  voltage  response,  where  statistical 
sampling  is  incorporated  to  account  for  the  system  uncertainties. 
This  yields  a  group  of  output  voltage  histories  as  the  baseline.  We 
finally  employ  the  Hotelling  T2  control  limit  approach  to  mon¬ 
itor  the  stack  condition  by  statistically  comparing  the  voltage 
measurement  with  the  baseline  voltages.  Our  analysis  indicates 
that  this  Hotelling  method  has  very  high  detection  sensitivity 
and  can  detect  the  fault  conditions  at  the  early  stage. 
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